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DEPENDENCE OF DUST OBSCURATION ON STAR FORMATION RATES IN GALAXIES 
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ABSTRACT 

Many investigations of star formation rates (SFRs) in galaxies have explored details of dust obscura- 
tion, with a number of recent analyses suggesting that obscuration appears to increase in systems with 
high rates of star formation. To date these analyses have been primarily based on nearby (z < 0.03) 
or UV selected samples. Using 1.4 GHz imaging and optical spectroscopic data from the Phoenix Deep 
Survey, the SFR-dependent obscuration is explored. The use of a radio selected sample shows that 
previous studies exploring SFR-dependent obscurations have been biased against obscured galaxies. The 
observed relation between obscuration and SFR is found to be unsuitable to be used as an obscuration 
measure for individual galaxies. Nevertheless, it is shown to be successful as a first order correction for 
large samples of galaxies where no other measure of obscuration is available, out to intermediate redshifts 
(zw0.8). 

Subject headings: galaxies: evolution — galaxies: starburst — radio continuum: galaxies 



1. INTRODUCTION 

Dust obscuration is currently recognized as one of the 
most serious sources of uncertainty in studies of galaxy 
evolution. With the recent results of far-infrared (FIR) 
and sub-mm observations revealing an ever increasing 
number of dusty star forming galaxies (e.g., Genzel & Ce- 
sarsky 2000; Ivison et al. 2000; Smail et al. 2002), the need 
to unify measures of star formation rate (SFR) from inde- 
pendent indicators at different wavelengths (e.g., UV, Ha, 
FIR, radio continuum) is as pressing as ever. 

A relatively simple prescription for dust extinction cor- 
rection to SFR has been suggested by Hopkins et al. (2001) 
and Sullivan et al. (2001), by assuming a luminosity- (or 
SFR) dependent obscuration. This was shown to provide a 
good first order correction to optically derived SFRs, while 
smaller differences still remain between different SFR in- 
dicators that are likely to be related to different star for- 
mation histories and/or extinction properties (Sullivan et 
al. 2001). 

The above studies were based on the comparison of dif- 
ferent SFR indicators for samples of relatively low redshift 
galaxies, selected at optical or UV wavelengths, which are 
prone to dust induced biases. Furthermore, at interme- 
diate redshifts where the Ha line falls out of the optical 
window and the [On]A3727 line can instead be used to 
measure the SFR (0.3 < z < 0.8), a direct comparison 
is more difficult (see Cardiel et al. 2003). In this paper 
we explore the validity of the luminosity- (or SFR) depen- 
dent obscuration at both low (z < 0.3) and intermediate 
(0.3 < z < 0.8) redshifts, using a sample of star-forming 
galaxies selected at radio wavelengths (which are insensi- 
tive to dust obscuration) with spectroscopic information. 
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2. OBSERVATIONS 



The Phoenix Deep Survey (PDS) includes a 1.4 GHz 
survey made using the Australia Telescope Compact Ar- 
ray (ATCA) covering a field a little over 4.5 square de- 
grees, selected to lie in a region of low optical obscura- 
tion and devoid of bright radio sources (Hopkins et al. 
2003a, and references therein). Optical imaging at the 
Anglo-Australian Telescope (AAT) has produced an op- 
tical catalogue probing to R = 22.5 for the whole field, 
and 2dF multi-object spectroscopy from the AAT provides 
spectra for many of the optically identified radio sources 
(Georgakakis et al. 1999; Afonso 2002). The spectra were 
taken using the low resolution gratings 270R, 316R and 
300B. Most objects were observed through only one of the 
red or blue gratings, with a small number being observed 
through both. The 270R and 316R gratings provide a 
wavelength coverage of 5000 < A < 8500 A, and the 300B 
gives 4000 < A < 7000 A. Redshifts were determined by 
visual inspection of the spectra, and line parameters were 
measured through Gaussian fitting using the SPLOT pack- 
age in IRAF. There are currently a total of 445 galaxies 
with measured spectra, of which 138 were securely classi- 
fied as star forming using optical emission line diagnostic 
diagrams. A more detailed account of the spectroscopic 
data reduction and classification can be found elsewhere 
(Georgakakis et al. 1999; Afonso 2002). 

3. STAR FORMATION RATES IN THE FAINT RADIO 
POPULATION 
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The origin of 1.4 GHz emission in star- forming galax- 
ies is primarily thought to be synchrotron radiation from 
relativistic electrons, accelerated by the shocks from su- 
pernova ejecta. The insensitivity of radio wavelengths to 
dust obscuration, makes radio emission a particularly at- 
tractive way of estimating SFRs in star-forming galaxies. 
The relation between SFR and 1.4 GHz luminosity can be 
written as 



SFRi .4 GHz 



Li. 



4 GHz 



8.4 x 10 20 W Hz 



— M yr" 



(1) 



for a Salpeter IMF with stellar masses between 0.1 — 
100 M Q (Haarsma et al. 2000). Using this relation, and 
calculating Li.4gh z assuming a radio spectral index a of 
0.8 (S„ oc v~ a ), characteristic of synchrotron radiation, 
the SFRs for the star-forming galaxies in the Phoenix 
spectroscopic sample were calculated and are presented 
in Figure 1. The depth and area of the radio survey and 
the spectroscopic follow-up allow sampling star-formation 
rates from below one to a few hundred solar masses per 
year. 

Of particular interest is the comparison between the 
SFRs as given by the radio and the optical line emission. 
To calculate the SFR from Ha, the conversion indicated 
by Kennicutt (1998), for an IMF and mass limits as for 
equation (1), is used 



SFRjja — 



1.27 x 10 34 W 



M yr- 1 . 



(2) 



The choice of a different IMF would change equa- 
tions (1) and (2) by a similar factor (which for the most 
commonly used IMFs would be ~2-4) , since both SFR in- 
dicators are sensitive to stars in roughly the same mass 
range (M>8M Q ). Hence, a comparison between the 
SFRs derived from radio and Ha emission will be insensi- 
tive to the particular IMF shape choosen. 

For galaxies at higher redshift (z > 0.3), where Ha 
could not be measured, one can use the [On]A3727 flux. 
This method relies on equation (2) and the observed ratio 
between [On] and Ha. A value of -Fjon] = 0.45 x Fn a 
(Kennicutt 1998) is commonly used. However, using the 
first data release of the Sloan Digital Sky Survey (Abaza- 
jian et al. 2003), Hopkins ct al. (2003b) observes a me- 
dian relation for radio-detected star-forming galaxies of 
Ffoni = 0.23 x Ffia (with a scatter of around 0.1), finding 
the difference to be due to the incompleteness of previ- 
ous samples. The range of optical luminosities sampled by 
Hopkins et al. (2003b) is similar to the one in the present 
Phoenix spectroscopic sample {Mr ~ 20 — 23.5) and is 
not enough to reveal the optical luminosity dependence in 
the [On] /Ha ratio observed by Jansen, Franx & Fabricant 
(2001), that produces lower [On] /Ha ratios for galaxies 
with higher optical luminosities. Also, the range of ex- 
tinctions, as given by the Balmer decrement, sampled by 
the Hopkins et al. (2003b) sample (Ha/H/3, between 3 and 
12) is similar to that of the present work (as will be seen 
below). We thus adopt the determination of Hopkins et 
al. (2003b) to convert our measured [On] luminosities to 
Ha values, using equation (2) to obtain the corresponding 
SFR estimate. 

Figure 2 compares the SFRs derived from the radio lu- 
minosity and line emission. Although a correlation exists, 



as expected, the SFRs derived from 1.4 GHz are in general 
higher that those calculated from nebular lines, especially 
for higher luminosities (or SFRs). This effect, seen pre- 
viously in several studies (Cram et al. 1998; Hopkins ct 
al. 2001; Sullivan et al. 2001), is attributed to dust obscu- 
ration, which affects the optical line emission. Further- 
more, one can conceive an amount of obscuration which 
increases with the SFR (a SFR-depcndcnt dust obscura- 
tion). This has been explored with considerable success 
for nearby optical (Hopkins et al. 2001) and UV selected 
galaxies (Sullivan et al. 2001). 



4. SFR-DEPENDENT DUST OBSCURATION 

For a subset of the Phoenix sample described above, 
an estimate of the extinction can be made from the ob- 
served Balmer decrement (Ha/H/3). Stellar absorption of 
the Balmer lines was corrected by assuming an average 
value of 2 A for the equivalent width (EW) of the H/3 ab- 
sorption in star-forming galaxies (Tresse et al. 1996; Geor- 
gakakis et al. 1999), with a similar value (2.1 A) being used 
for the Ha line (equation (2) of Miller & Owen 2002). Fig- 
ure 3 shows the resulting Balmer decrements, corrected for 
stellar absorption, as a function of SFR, derived from the 
radio luminosity using equation (1). Unlike previous stud- 
ies (Hopkins et al. 2001; Sullivan et al. 2001), a tight corre- 
lation is not observed here. Rather, a trend for a broader 
range of obscurations for higher SFR systems is observed. 
This behaviour seems to be due to different selection cri- 
teria and small number statistics, as we now explain. 

The sample used by Hopkins ct al. (2001) to derive 
the relation between SFR and obscuration includes only 
nearby (z < 0.03) galaxies with EW(Ha) larger than 30 A. 
While no clear trend is seen in the present sample when 
restricted to this EW value, the higher limit of 60 A (filled 
circles in Figure 3) does suggest a closer match to the ob- 
servations in Hopkins et al. (2001). 

On the other hand, Buat et al. (2002) also observed a 
dual behaviour: while nearby star forming galaxies behave 
similarly to what is present in Figure 3, a sample of IUE 
galaxies shows a much tighter correlation, as that observed 
for the UV-selected sample of Sullivan et al. (2001). This 
suggests that UV-selection results in some kind of bias 
that is avoided with the present sample. Completeness in 
Figure 3 is not easy to quantify, given the several selection 
criteria (initially the radio flux limit, followed by the op- 
tical identification and 2dF spectroscopy, which imposes 
a practical limit of R <~ 20). However, it is possible to 
try to understand a possible bias in a magnitude-limited 
UV selected sample, and at the same time to evaluate the 
improvement of the present work, as we now show. 

The tight relation between dust-free UV emission and 
SFR can be used to evaluate which regions of Figure 3 
are not accessible to a magnitude-limited UV study. As- 
suming a limiting magnitude tojjv — 18-5 (as in Sullivan 
et al. 2001), an intrinsic SFR (ic, L\jy before obscura- 
tion) and redshift will define the maximum value of the 
Balmer decrement that still allows a detection. Figure 4 
shows the present sample, separated according to redshift, 
overlaid with the maximum detectable Balmer decrement 
at z = 0.05 (dotted line), 0.1 (dashed line) and 0.2 (dot- 
dashed line). The conversion between SFR and L\jy uses 
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the calibration from Kennicutt (1998), while the extinc- 
tion at 2000 A is derived from the Balmcr decrement us- 
ing the procedures of Calzetti et al. (2000). An estimate 
of the X-correction is obtained using an average colour of 
m uv - b = -1.5 (Milliard et al. 1992). It is clear that 
the present sample represents a significant improvement 
for z > 0.1. In particular, many of the galaxies showing 
high Balmer decrement values in the present study would 
not be detected in a UV survey limited to myv = 18.5. 
Sample selection thus seem to be a major source of bias 
when trying to investigate the correlation between dust 
obscuration and SFRs. 

Given the large scatter present in Figure 3, a SFR- 
dependcnt reddening correction is obviously unsuitable for 
application in galaxies where a direct estimate of obscu- 
ration exists. However, a trend for higher average Balmer 
decrement (and greater distribution width) with increasing 
SFR seems to exist. This can still be useful as a prelimi- 
nary dust obscuration estimate for large samples of galax- 
ies where no other measure of obscuration is available. Al- 
though in practice the form of the derived relation may be 
comparable to the ones in Hopkins et al. (2001) and Sulli- 
van et al. (2001), here we recognize that there is no tight 
correlation between obscuration and SFR, but an average 
obscuration may still be defined for any given SFR. As can 
be seen in Figure 3 the resulting correction will be affected 
by large uncertainties for individual galaxies, especially at 
large SFRs. 

The sample was thus split into 7 bins of log(SFR) (as es- 
timated from the radio luminosity) , each having between 5 
and 16 objects. The median log(SFR) and Balmer decre- 
ment in each bin were then found (shown as asterisks in 
Figure 3). A linear fit, taking into account the errors in 
both quantities, results in 



Ha 



77^) = 1.29 log(SFR) + 5.06, 

i / median 



(3) 



with a correlation coefficient of 0.8. Keeping in mind the 
meaning and limitations of this correlation, as seen in Fig- 
ure 3, one can now test its usefulness as a first correction 
for the effect seen in Figure 2. 

The departure of the observed Balmer decrement from 
the Case B value of 2.86 (e.g., Brocklchurst 1971), can 
be related to the color excess for nebular emission lines, 
E(B — V) gas , and extinction, k(X), by 



HtA =(^£L\ -, () 0AE(B-V) ga Ak a 



(4) 



obs 



Substituting (3) into (4) gives a relation for the color ex- 
cess as a function of SFR: 



E(B - V) gas = 



2.5 



log 



{k a -k ) b Vl.291og(SFR) + 5.06 



(5) 

Together with an appropriate extinction curve (the stan- 
dard Galactic extinction curve of Cardclli, Clayton & 
Mathis (1989) with R v = 3.1, found by Calzetti (2001) 
to describe well the reddening of the ionized gas in star- 
forming galaxies), this can then be used to correct Lu a , 
and consequently, SFRh q , for dust obscuration: 



Lua = it 10 4 E{B - k " . (6) 

where Lg^ can either be the observed Ha luminosity or 
the "effective" Ha luminosity derived from an observed 
[On] flux. 

Equation (5) gives the relation between extinction and 
the intrinsic SFR. Assuming this to be the value given 
by the radio luminosity could be a good approximation, 
but would create an artificial dependence between the cor- 
rected Ha SFR and the one from 1.4 GHz. Instead, since 
the form for the SFR-dependent obscuration is monoton- 
ically increasing, an iteration over possible values for in- 
trinsic SFR and the corresponding obscuration can be per- 
formed until the calculated obscured SFR converges with 
the observed value (Hopkins et al. 2001). We note that 
this procedure does not take into account any absorption 
of ionizing photons by dust inside HII regions. Chariot et 
al. (2002), modeling the observed spectra in non-Seyfert 
galaxies, estimate that this mechanism is responsible for 
the loss of ^20% of ionizing photons. Given the large un- 
certainty associated with this value, however, we do not 
attempt any correction, noting that its magnitude would 
not significantly affect our results. 

Figure 5 shows the resulting dust corrected relation for 
the SFR from line and radio luminosities. It is clear that 
the SFR-dependent dust absorption, while being a very 
coarse approximation, can successfully account for the first 
order offset between the SFRs derived from Ha or [On] and 
radio luminosities for galaxies spanning a broad range of 
rcdshifts (out to z w 0.8). This would not be possible if 
the relations between Balmer decrements and SFR drawn 
from previous samples (Hopkins et al. 2001; Sullivan et 
al. 2001) had been used. The scatter still present has an 
rms of 0.4 dex about the best fit line, maintained from the 
scatter in Figure 2. The lack of an improvement lies in 
the coarse relationship between SFR and obscuration seen 
in Figure 3 - the linear fit to the median values cannot 
correct for the range of obscurations seen at each SFR. 

There will be, of course, additional uncorrelated mech- 
anisms involved in the Ha and radio emission which con- 
tribute to the scatter seen, but their quantification will 
only be possible after a precise account of the obscuration 
for each individual galaxy. 

5. SUMMARY 

A radio selected sample of star forming galaxies to 
z ss 0.8 has been compiled from the Phoenix Deep Survey. 
The use of radio selection minimises bias in the sample 
due to dust obscuration effects. The relationship between 
obscuration and SFR is shown to be only for a higher 
Balmer decrement range at higher SFRs, contrary to the 
tight correlation observed in previous studies. Still, the 
use of a linear relation, reflecting only the broadest trend, 
was explored as a first order correction for large samples of 
galaxies with no direct measurement of obscuration. This 
successfully accounts for the major discrepancy between 
optical emission line SFR estimates and 1.4 GHz luminos- 
ity estimates for all galaxies in the broad redshift range 
probed. However, a much more detailed correction of the 
dust obscuration is necessary for the study of the uncorre- 
lated mechanisms (e.g., star formation histories) responsi- 
ble for the scatter still present. 
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Fig. 1. — SFRs for star-forming galaxies in the Phoenix Deep Survey, as estimated from their radio luminosity. Crosses indicate galaxies 
with an Ha measurement suitable for a SFR estimate, while circles denote those where only [Oll]A3727 is measured. The dashed line indicates 
the 60 /xJy detection limit of the radio observations. 
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Fig. 2. — Comparison between the SFRs derived from 1.4 GHz and from Ha or [On] luminosities. The crosses represent galaxies with direct 
Ha measurement, circles denote galaxies at higher redshift, where [On] was used as a proxy for the Ha line intensity. 
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Fig. 3. — Balmer decrements, corrected for stellar absorption, as a function of SFR derived from radio luminosities. The solid line indicates 
the best fit to the median values (asterisks) in bins of log(SFR), while the dotted line shows the relationship derived in Hopkins et al. (2001) 
for a sample of local galaxies. The error bar in the bottom left corner shows the median errors for the Balmer decrement and radio flux 
measurements. 



8 



Afonso et al. 



15 



i — i — i i 1 1 



T 1 1 1 — I I I I 



- 1/ 1 1 1 — I I I'lj 

/z=0.05 /z=0.1 



*z < 0.05 
o 0.05 < z ^ 0.1 
• 0.1 < z g 0.2 
*z > 0.2 



10 



/ 

/ 

• / * 

• */ 

/ 

/ • 



■P.-' 



/ 

* * 
/ * 



/ 

* / 

/ 

/' 

/' 

/ 

/z=0.2 J 



o 



/ 

/O 



t 



.'it 



/ 



it 



/ 



Q» y'&t o*„ 



05 



J 



J I I 



J I I 



10 

SFR 1 . 4GHz (M yr-0 



100 



Fig. 4. — Data as in Fig. 3, marked according to redshift. The lines indicate the maximum Balmer decrement detectable for a given SFR 
at redshifts of 0.05, 0.1 and 0.2, in a UV survey with a limiting magnitude of rrtyv = 18.5. It is clear that most of the galaxies detected here 
displaying the highest Balmer decrements, would not be detected in such a UV survey. 
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. 5. — Same as Fig. 2, but with line luminosities corrected for dust obscuration using Equation (6). 



